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Abstract 


Classical  Ekman  spiral  is  generated  by  surface  wind  stress  with  constant  eddy  viscosity  in 
homogeneous  ocean.  In  real  oceans,  the  eddy  viscosity  varies  due  to  turbulent  mixing  caused  by 
surface  wind  and  buoyancy  forcing.  Horizontally  inhomogeneous  density  produces  vertical 
geostrophic  shear  which  contributes  to  current  shear  that  also  affects  the  Ekman  spiral.  Based  on 
the  similar  theoretical  framework  as  the  classical  Ekman  spiral,  the  baroclinic  components  of  the 
Ekman  spiral  caused  by  the  horizontally  inhomogeneous  density  are  obtained  analytically  with 
the  varying  eddy  viscosity  calculated  from  surface  wind  and  buoyancy  forcing  using  the  K- 
Profile  Parameterization  (KPP).  Along  with  the  three  existing  types  of  eddy  viscosity  due  to 
pure  wind  forcing  (zero  surface  buoyancy  flux),  such  an  effect  is  evaluated  using  the 
climatological  monthly  mean  data  of  surface  wind  stress,  buoyancy  flux,  ocean  temperature  and 
salinity,  and  mixed  layer  depth. 
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1.  Introduction 


On  the  base  of  homogeneous  density  without  considering  waves,  Ekman  (1905)  modeled 
turbulent  mixing  in  upper  ocean  as  a  diffusion  process  similar  to  molecular  diffusion,  with  an 

eddy  viscosity  (turbulent  plus  molecular),^  (the  symbol  6A’  indicating  dimensional  quantity), 
which  was  taken  as  a  constant  with  many  orders  of  magnitude  larger  than  the  molecular 
viscosity.  The  turbulent  mixing  generates  ageostrophic  component  of  the  upper  ocean  currents 
(called  the  Ekman  spiral),  decaying  by  an  e-folding  over  a  depth  as  the  current  vector  rotate  to 
the  right  (left)  in  the  northern  (southern)  hemisphere  through  one  radian.  Several  approaches  may 
advance  the  classical  Ekman  theory:  (a)  replacing  constant  eddy  viscosity  by  varying  eddy 
viscosity,  and  relating  the  eddy  viscosity  to  ocean  mixing  (under  surface  wind  and/or  buoyancy 
forcing),  (b)  including  ocean  wave  effect,  and  (c)  changing  homogeneous  to  inhomogeneous 
density. 

It  was  recognized  that  the  eddy  viscosity  K  is  not  a  constant.  After  fitting  observational 
ocean  currents  to  the  Ekman  spiral  (e.g.,  Hunkins  1966;  Stacey  et  al.  1986;  Price  et  al.  1987; 

Richman  et  al.  1987;  Chereskin  1995;  Lenn  and  Chereskin  2009),  the  inferred  K  value  varies 
more  than  an  order  of  magnitude,  from  0.054  m  s'  (Price  et  al.  1987)  obtained  from  the  field 
measurements  acquired  from  a  surface  mooring  set  in  the  western  Sargasso  Sea  (34°N,  70°W)  as 
part  of  the  Long  Term  Upper  Ocean  Study  Phase  3  (LOTUS-3)  during  the  summer  of  1982,  to 
0.006  m  s'  (Stacey  et  al.  1986)  obtained  from  the  low-frequency  current  measurements  in  the 
Strait  of  Georgia,  British  Columbia.  The  smaller  value  (0.006  ms')  may  be  treated  as  a  lower 
bound  of  the  eddy  viscosity  (Price  et  al.  1987).  Recently,  Lenn  and  Chereskin  (2009)  obtained  a 
mean  Ekman  spiral  from  high-resolution  repeat  observations  of  upper-ocean  velocity  in  Drake 
Passage  along  with  the  constant  temperature  in  the  Ekman  layer  (implying  near  neutral 
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stratification).  The  eddy  viscosities  inferred  from  Ekman  theory  and  the  time-averaged  stress  was 
directly  estimated  as  O(10"2-10_1)  m2  s"1. 

The  turbulent  mixing  in  upper-ocean  is  also  viewed  as  being  driven  by  the  atmospheric 
fluxes  of  momentum  and  buoyancy  (heat  and  moisture),  and  the  shear  imposed  by  the  ocean 
circulation,  and  characterized  by  the  existence  of  a  vertically  quasi-uniform  layer  of  temperature 
and  density  (i.e.,  mixed  layer).  Underneath  the  mixed  layer,  there  exists  another  layer  with  a 
strong  vertical  gradient,  such  as  the  thermocline  (in  temperature)  and  pycnocline  (in  density) 
(e.g.,  Kraus  and  Turner  1967;  Garwood  1977;  Chu  and  Garwood  1991;  Steger  et  al.  1998;  Chu  et 
al.  2002).  Such  vertical  mixing  generates  varying  upper-ocean  eddy  viscosity.  The  mixed  layer  is 
a  key  component  in  studies  of  climate  and  the  link  between  the  atmosphere  and  deep-ocean  and 
directly  affects  the  air-sea  exchange  of  heat,  momentum,  and  moisture  (Chu  1993). 

Effect  of  vertical  inhomogeneity  of  density  on  the  Ekman  spiral  (i.e.,  stratified  Ekman 
layers)  has  been  identified  by  observational  and  modeling  studies  in  atmospheric  boundary  layer 
(Lettau  and  Dabberdt  1970;  Grachev  et  al.  2008)  and  oceanic  boundary  layer  (McWilliams  et  al. 
2009;  Taylor  and  Sarkar  2008).  Ocean  observations  from  drifter/floats  show  the  role  of 
horizontal  density  gradient  in  setting  the  stratification  within  the  mixed  layer.  McWilliams  et  al. 
(2009)  computed  vertical  turbulent  mixing  within  the  boundary  layer  in  a  one-dimensional 
vertical  column  using  the  K-profile  parameterization  (KPP)  scheme  with  surface  mean  wind 
stress,  mean  heating,  and  solar  absorption,  and  idealized  representations  of  the  heat  flux  from  the 
interior  three  dimensional  circulation  and  found  that  there  is  not  a  single,  simple  paradigm  for  the 
upper-ocean  velocity  profiles  in  stratified  Ekman  layers  due  to  the  following  reasons:  (a)  the 
Ekman  layer  is  compressed  by  stable  stratification  and  surface  heating;  (b)  Ekman  currents 
penetrate  down  into  the  stratified  layer;  (c)  penetrative  solar  absorption  deepens  the  mean  Ekman 
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layer;  (d)  wind  and  especially  buoyancy  rectification  effects  yield  a  mean  Ekman  profile  with  a 
varying  eddy  viscosity,  where  the  mean  turbulent  stress  and  mean  shear  are  not  aligned,  whereas 
buoyancy  rectification  induces  profile  flattening.  These  modeling  results  are  for  the  one¬ 
dimensional  ocean,  i.e.,  no  horizontal  gradients  of  any  variables  including  the  density. 

Effect  of  ocean  surface  gravity  waves  on  the  Ekman  spiral  has  been  identified  through 
interacting  waves  with  ocean  currents  and  wind  stresses.  As  waves  experience  breaking  and 
dissipation,  momentum  passes  from  waves  into  ocean  currents.  Recent  studies  show  that  the 
influence  of  the  surface  wave  motion  via  the  Stokes  drift  and  mixing  is  important  to 
understanding  the  observed  Ekman  current  profiles  in  addition  to  wind  stress,  depth-varying 
eddy  viscosity,  and  density  inhomogeneity.  Song  and  Huang  (2011)  used  the  WKB  method  to 
obtain  the  analytic  solutions  for  modified  Ekman  equations  including  random  surface  wave 
effects  when  the  eddy  viscosity  is  gradually  varying  with  depth.  Their  solution  was  compared 
with  observational  data  and  with  the  results  from  a  large  eddy  simulation  of  the  Ekman  layer 
(Zikanov  et  al.  2003). 

However,  effect  of  horizontally  inhomogeneous  density  on  the  Ekman  spiral  with  varying 
eddy  viscosity  due  to  vertical  mixing  under  various  surface  forcing  conditions  has  not  yet  been 
studied.  Since  horizontally  inhomogeneous  density  leads  to  non-zero  vertical  geostrophic  shear, 
and  in  turn  contributes  to  the  current  shear,  the  equations  and  surface  boundary  conditions  for  the 
classical  Ekman  model  need  to  be  modified.  Such  modifications  may  lead  to  a  new  structure  of 
the  Ekman  spiral.  The  baroclinic  components  of  the  Ekman  spiral  are  identified  analytically  in 
this  study  using  the  KPP  and  three  existing  (due  to  pure  wind  forcing)  eddy  viscosities  without 
considering  ocean  waves.  Rest  of  the  paper  is  organized  as  follows.  Section  2  introduces  the 
basic  equations  and  boundary  conditions.  Sections  3  and  4  describe  the  Obukhov  length  scale 
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(Obukhov  1946,  Monin  and  Obukhov  1954),  depth  ratio,  and  KPP.  Section  5  presents  the 
analytical  solution  of  the  Ekman  spiral  in  horizontally  inhomogeneous  ocean  including  analytical 
barotropic  and  baroclinic  components  due  to  KPP  eddy  viscosity.  Section  6  and  7  describe  the 
baroclinic  effects  with  the  KPP  eddy  viscosity  under  both  surface  wind  and  buoyancy  forcing 
and  with  the  three  existing  eddy  viscosities  under  pure  wind  forcing.  Section  8  presents  the 
conclusions.  Appendices  A  and  B  list  the  procedures  for  obtaining  the  analytical  solutions  of  the 
Ekman  spiral  in  horizontally  inhomogeneous  ocean  with  depth-dependent  eddy  viscosity. 

2.  Ekman  Layer  Dynamics 


Let  (x,  y9  z )  be  the  zonal  (positive  eastward),  latitudinal  (positive  northward),  and  vertical 
(positive  upward  with  z  =  0  at  the  ocean  surface)  coordinates  with  (i,  j,  k)  as  the  corresponding 
unit  vectors,  and  u  be  the  velocity  vector.  Following  the  similar  steady  dynamics  of  McWilliams 
and  Huckle  (2006)  with  modification  from  homogeneous  to  inhomogeneous  density,  the  steady- 
state  horizontal  momentum  balance  with  Boussinesq  approximation  is  given  by 


/|  /V  1  „  is 

ykxu  = - Vp— -  — 

Pw  hdcr 


(la) 


where  pw  =  1025  kg  m-3,  is  the  characteristic  density  of  seawater;  h  is  the  ocean  surface  mixed 

layer  depth;  a  =  -z  /  h ,  is  the  non-dimensional  vertical  coordinate;  /  is  the  Coriolis  parameter 
(depending  on  the  latitude);  M  is  the  vertical  momentum  flux  due  to  turbulent  mixing;  p  is  the 
pressure.  It  is  noted  that  the  damping  for  currents  due  to  vertical  radiation  of  inertial  waves  into 
the  oceanic  interior  is  neglected.  The  mixed  layer  depth  ( h )  can  be  determined  from  temperature 
and  density  profiles  using  subjective  and  objective  methods  (e.g.,  Monterey  and  Levitus  1997; 
Chu  et  al.  2002;  Chu  and  Fan  2010,  2011).  The  hydrographic  balance  gives 


1  dp  p 

—  —  =  g~g  — , 

A,  &  pw 


(lb) 
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where  p  is  the  density;  g  is  the  gravitational  acceleration  (9.81  m  s"2). 

The  horizontal  velocity  consists  of  two  parts:  geostrophic  current,  [u  =  (^g,vg)],  and 


ageostrophic  current  [u^  =  (uE,vE)  ]  (Ekman  flow), 


u  =  ug+uE. 


(2) 


where  the  geostrophic  current  is  given  by 


1 


ykxu  = — vp, 

pw 


(3) 


and  computed  solely  from  the  density  field  (Chu  1995;  2000;  2006).  Differentiation  of  (3)  with 
respect  to  z  and  use  of  (lb)  lead  to  the  thermal  wind  relation 

5u„ 


dz 


=  -kx 


g 


■Vp 

\fPw  j 


(4) 


Substitution  of  (2)-(3)  into  (la)  leads  to 


»  ~  15 

/kxu,  = - 

£  hda 


The  vertical  momentum  flux  M  (i.e.,  turbulent  Reynolds  stress)  is  modeled  by 


x  K  5u 
M(z)  =  — , 
h  da 


where  K  is  the  eddy  viscosity  that  is  non-dimensionalized  by 


(5) 


(6a) 


K  = 


K 

h/cu * 


(6b) 


where  k~0.4  1 ,  is  the  von  Karmen  constant.  Substitution  of  (2)  into  (6a)  and  use  of  (6b)  lead  to 


5u  (3u 

M(z)  =  -ku,K  ^  +  ^ 
dcr  da 


(V) 


The  velocity  (uE )  is  nondimesionalized  by 
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(8) 


Ue  = 


vK 


VP  = 


where  7r(0)  is  the  eddy  viscosity  evaluated  at  the  surface.  Substitution  of  (7)  into  (5)  and  use  of 


(4)  and  (8)  give 


JkxuE- 


kiu.  d 


h  da 


K ^ 

da 


tcu*  ,  d 


-kx  —  \KS],  S  =  -^-Vp 

JVE  J  pw 


(9) 


Here,  the  vector,  S  =  (sx,  sy),  is  defined  by 

(10) 

Pw  dx  pw  dy 

which  represents  the  baroclinicity  (i.e.,  sx  ^  0,  s  *  0 ).  The  ocean  is  barotropic  if 

5X=  Sy  =  0.  (11) 

The  second-order  differential  equation  (9)  needs  two  boundary  conditions.  At  the  surface 
(cr  =  0)  we  have 


M(0)  =  —  =  CD  |ufl | ua  =  w»0» 

Pw  Pw 


(12a) 


/V  ^ 

where  Q>  is  the  drag  coefficient;  t  is  the  surface  wind  stress;  pa  =  1.29  kg  m'  ,  is  the 
characteristic  atmospheric  density;  ua  is  the  wind  near  the  ocean  surface;  0,  =  [cos  6,  sin  6\  is 

the  unit  vector  of  the  wind  direction;  6  is  the  angle  of  the  wind  from  the  east;  and  u*  is  the 
ocean  friction  velocity, 


U  A 


^DPa  |Uq| 

Pw 


1/2 


(12b) 


Evaluation  of  (7)  and  (4)  at  the  surface  leads  to 


M(0)  =  -ku„K 


^ 8  +V»e 


da  da 


(13a) 


y  cr=0 
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(13b) 


f 

=  kx 

(  7  A 

g  Vp 

cr=0 

Pw  y 

Substitution  of  (10),  (12a),  and  (13b)  into  (13a)  leads  to  the  surface  boundary  condition  for  the 
non-dimensional  Ekman  flow  u^, 


V  dcF  Ja=Q 


u * 


0*  - 


h 


kK(0)Ve  JVe 


kxS(0). 


(14) 


finite  as  cr^oo 


where  K(0)  and  S(0)  represent  the  values  of  (K,  S)  evaluated  at  the  surface  (o  =  0).  Moreover,  the 
general  solution  of  (9)  contains  exponentially  increasing  and  decreasing  parts  with  the  non- 
dimensional  depth  a.  The  exponentially  increasing  part  is  unphysical  and  needs  to  be  eliminated. 
Therefore,  the  lower  boundary  condition  of  equation  (9)  is  used 

(15) 

to  filter  out  the  unphysical  solution.  In  fact,  the  lower  boundary  condition  (15)  is  also  used  in  the 
classical  Ekman  spiral. 

Generally,  equation  (9)  is  not  closed.  One  more  equation  for  the  density  p  is  needed.  If  p  is 
given,  the  second-order  differential  equation  (9)  with  the  boundary  conditions  (14)  and  (15)  are 
well-posed.  For  depth-dependent  eddy  viscosity,  (9)  is  an  inhomogeneous  linear  differential 
equation  with  variable  coefficient  K. 

3.  Obukhov  Length  ( L )  and  Depth  Ratio  (>.) 

The  eddy  viscosity  is  to  characterize  vertical  mixing,  which  is  generated  by  surface  wind 
stress  (t)  and  surface  buoyancy  flux  ( B  in  m2  s"3,  upward  positive), 


B  =  -^=-  +  g(3{E-P)S,  (16) 

P(0)cp 

where  Q  is  the  net  heat  flux  (upward  positive,  W  m' );  cp  is  the  specific  heat  for  the  sea  water;  S 
is  the  surface  salinity  [in  practical  salinity  units  (psu)];  a  is  the  coefficient  of  thermal  expansion; 
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P  is  the  coefficient  of  haline  contraction;  and  (E,  P )  are  evaporation  and  precipitation  (m  s"1). 
Ocean  mixed  layer  is  generally  developed  by  wind  stirring  and  convection  (upward  surface 
buoyancy  flux  B ).  To  examine  dominant  mixing  mechanisms,  the  Obukhov  length  scale  ( L )  and 
the  depth  ratio  (A)  are  calculated  by 

Z=-iA  /t=*=a0m)  h  .  (17) 

kB  L  L  (10  m) 

Here,  L  is  the  depth  where  the  wind-generated  turbulence  is  balanced  by  the  downward 
buoyancy  flux  (B  <  0)  due  to  surface  warming  ( Q  <  0)  and/or  freshening  (P  >  E)  and  is 
comparable  to  the  convection-generated  turbulence  by  the  upward  buoyancy  flux  ( B  >  0)  due  to 
surface  cooling  ( Q  >  0)  and/or  salinization  (P  <  E);  and  A  is  the  depth  ratio.  Monthly  depth  ratio 
(X)  (Fig.  1)  are  calculated  from  the  monthly  mean  global  ocean  (10  m )/L  and  ocean  friction 
velocity  (u*)  data  (1°  x  1°)  downloaded  from 

http://iridl.ldeo.columbia.edU/SOURCES/.DASILVA/.SMD94/.climatology/  (DaSilva  et  al. 
1994),  and  the  monthly  mixed  layer  depth  (h)  data  downloaded  from 
http://www.nodc.noaa.gov/OC5/WOA94/mix.html  (Monterey  and  Levitus  1997). 

The  depth  ratio  (A,)  is  used  to  determine  the  forcing  regimes  (Lombardo  and  Gregg  1989): 
convective  regime  ( A<  -10),  wind-forcing  regime  (A>- 1),  and  combined  forcing  regime 
(-10  <  A  <  —  1 ).  The  depth  ratio  (A)  also  serves  as  a  stability  parameter  (see  next  section).  The 
calculated  monthly  depth  ratio  (A)  (Fig.  1)  shows  strong  seasonal  variability  with  only  two 
regimes  evident:  wind- forcing  and  combined  forcing  regimes  since  almost  no  data  with  A  <  -10. 
In  January,  the  combined  forcing  ( — 10  <  A  <  — 1 )  prevails  most  of  the  northern  hemisphere 
including  North  Atlantic,  North  Pacific,  Arabian  Sea,  Mediterranean  Sea,  and  eastern  tropical 
South  Pacific;  and  the  wind  forcing  (A  >  -1 )  prevails  most  of  the  southern  hemisphere.  In  July, 


10 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 


the  combined  forcing  prevails  in  the  southern  hemisphere;  and  the  wind  forcing  prevails  in  the 
northern  hemisphere. 

4.  KPP 


With  the  surface  wind  and  buoyancy  forcing,  the  KPP  rules  for  the  non-dimensional  eddy 
viscosity  ( K)  are  given  by  the  product  of  a  depth-dependent  non-dimensional  turbulent  velocity 
wx(o)  (scaled  by  Km)  and  a  dimensionless  vertical  shape  function  G{o)  (Large  et  al.  1994) 


K(cr,A)  = 


fwx(cr,/l)G(cr),  if  1  >  cr  >  0 


wx(U)G(  1), 


if  <j>\  ’ 


(18) 


to  represent  the  capability  of  deeper  mixed  layers  to  contain  larger  more  effect  turbulent  eddies. 
It  is  noted  that  the  extension  of  eddy  viscosity  K(a,  A)  from  1  >  cr  >  0  to  <x  >  0  is  due  to  the  fact 
that  h  was  defined  as  the  boundary  layer  depth  in  the  original  KPP  model,  which  is  usually 
deeper  than  the  mixed  layer  depth.  The  shape  function  G(a)  is  assumed  to  be  a  cubic  polynomial 
(O’Brien  1970)  and  given  by  (McWilliams  and  Huckle  2006)  (Fig.  2) 


G(cr)  =  cr  [l  -  cr2  ]  +  ^  —  H(a0  -  cr), 

L  J  2crn 


(19) 


where  <x0  =  0.05,  H(a)  is  the  Heaviside  step  function  (equal  to  1  for  a  >  0  and  0  otherwise).  As 

pointed  out  by  McWilliams  and  Huckle  (2006),  the  second  term  in  the  righthand  side  of  (19)  is 
the  mathematical  aesthetics  and  computational  regularity. 

The  depth  dependent  non-dimensional  turbulent  velocity  scale  wx(a)  is  given  by  (Large  et 
al.  1994) 


wx(o-,A)  =  < 


1 


<p{sX) 


<j){crX) 


,  £<cr<l  T<0,  £=0.1 


otherwise 


(20) 
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Here,  the  function  </>  is  defined  by  the  Monin-Obukhov  similarity  theory  (Monin  and  Obukhov 
1954)  such  that  the  dimensional  turbulent  velocity  scales  equal  ku*  with  neutral  forcing  (k  =  0) 
and  are  enhanced  and  reduced  in  unstable  (k  <  0)  and  stable  ( k  >  0)  conditions.  It  is  given  by 
(Large  et  al.  1994) 


(/>(crA)  = 


1  +  5a>l,  0  <  A 

(l-16aAy1/4,  A<0,  cr<am/A 
(1.26-8.38a>l)“1/3,  A<0,  a  >am/  A 


am  =  -0.2 


(21) 


For  neutral  forcing  (X  =  0),  <f>(crA)  =  1 .  Substitution  of  (19),  (20),  and  (21)  into  (18)  leads  to  an 
analytical  non-dimensional  KPP  eddy  viscosity  K(a,  A).  For  a  given  depth  a,  <f>  increases  with  A 
(Fig.  3a);  and  K(a,  X)  decreases  with  A  (Fig.  3b).  Such  2-dependence  of  (j)  and  K(a,  X)  is  quite 
smooth  for  A  >  0  and  A  <  0,  but  very  abrupt  at  2  =  0.  The  <f)  -  values  are  small  for  A  <  0  (e.g., 


(j)  =  0.05  for  a  =  0.5,  A  =  -1)  and  very  large  for  A  >  0  (e.g.,  (j)  ~  30,000  for  a  =  0.5,  A  =1).  The  K- 


values  are  large  for  A  <  0  (e.g.,  K  =  7  for  a  =  0.5,  A  =  -1)  and  very  small  for  A  >  0  (e.g.,  K  « 10  5 
for  a  =  0.5,  A  =1).  However,  the  dependence  of  (/)  and  K(a,  X)  on  a  is  quite  mild.  Substitution 
of  the  KPP  eddy  viscosity  at  the  surface  K( 0,  A)  into  (8)  leads  to 

uyl 

V  = _ _ _ 

p\f\hKK(0,A) 

Monthly  Ekman  velocity  scale  ( VE),  calculated  from  the  same  data  sets  for  the 
computation  of  2,  has  strong  seasonal  variability  (Fig.  4).  In  January,  larger  f^-values  (>  0.5  m 
s'1)  occur  in  the  northern  hemisphere  such  as  in  the  Gulf  Stream,  Kuroshio,  equatorial  regions 
(especially  in  the  eastern  Pacific),  and  smaller  f^-values  (<  0.2  m  s'1)  occur  in  the  southern 
hemisphere.  In  July,  smaller  F^-values  (<  0.2  m  s'1)  occur  in  the  northern  hemisphere  except 
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some  northern  tropical  regions  such  as  near  the  northern  African  coast  and  west  Arabian  Sea,  and 
larger  F^-values  (>  0.5  m  s"1)  occur  in  the  southern  hemisphere. 

5.  Ekman  Spiral 


Substitution  of  (18)  into  (9)  leads  to 


ykxu  *- 


ku*  d 


h  da 


K(c 7,  A)-5"* 


KU*- k  x  [K(a ,  /l)S] , 


d<y  J  jVE  dc 


(22) 


which  is  an  ordinary  differential  equation  with  depth-varying  K(o,  X).  The  WKB  method  was 
used  in  this  study  to  solve  the  differential  equation  (22)  with  the  boundary  conditions  (14)  and 
(15)  to  get  the  approximate  analytical  solutions  uE  [=  (uE,  v/,;)]  (see  Appendix  A), 


Mr  Wr  A  Mr ,  vE  Vp  -I-  Avy , 


(23) 


where 

uE  =  exp  [F(cr,  /I)]  |a+  cos  [F(cr,  /l)]  +  a~  sin  [F(cr,  /l)]  j ,  (24a) 

vE  =  exp \F(cr,  /I)]  | a+  sin  [F(cr,  A)]  -  a~  cos  [F(cr,  /l)]}  (24b) 


are  the  barotropic  components  of  the  Ekman  velocity  (i.e.,  sx  =  0 ,sy  =  0);  and 
A uE  =  y  exp  [F(< j,  A)]  j  J3+  cos  [F(cr,  A)\  +  /T  sin  [F(cr,  T)]| 


Y  sgn(/) 

2  f2  ' 


8[4k<YT)(s,  -j,)] 
&T 


x  y  dC 


cos  [F(cj,A)-F(<Z,A)\ 


+  (^+^) 


ay  k(£,a) 
ac 


sin  [F(*,A)-F(£,A)] 


c  exp  ([^X<r ,  A)  -  F(£ ,  A)]) 
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rsgn(/) 

2/2  ; 


4mJ )(*,-*,)' 


+ 


T) 

*  y  dc 


cos \F(c t,  A)  -  F(£,  A)] 


dC 


+  6*+^) 


sm[F(a,A)-F{<;,Aj\ 


-  exp  (- [F(cr,  A)  -  F(£,  A)\) 


Ave  =  y  exp  [F(<t,  Aj\  j/?+  sin  [F(cr,  A)]  -  fd  cos \F(cj,  A)]| 

djK(£,A) 

df 

djK (£,A) 


rsgn(/) 

2  f2 


c 

i  [4mJ)(Sx+sr) 

d£ 

V 

A 

i 

d[4mJ)(S,-Sy) 

U 

+  (SX+Sy) 


A  +  (sx-sv) 


dc 


cos  [F(ct,A)-F(<Z,A)] 


sin  [F(cr ,  A)  -  F{C, ,  A)] 


c  exp  ([F(cr,  A)  -  F{£,  A)]) 


Y  sgn(/) 

2 /2  : 


\m^)(sx+Sy) 


dC 

7  W^)(sx-sy) 


+  (SX+Sy) 


+  (Sx~Sy) 


djKj £,A) 

djKA^F) 


cos \F(ct,  A)  -  F(£,  A)] 


sin  [T7  (cr,A)-  F(£ ,  A)\ 


:  exp  (-  [F(c t,  A)  -  F(£,  Aj\) 


are  the  baroclinic  components  of  the  Ekman  velocity  (i.e.,  nonzero  if  sx  *  0,  s 
parameters  are  defined  as  follows: 


Kh\f\JK(0,X)  +  fl,  if  /  >  0 

7  = — — - ,  a~  =(cos6»±sin6»)sgn(/),  sgn(/)  =  j 

M*  [-1,  if  /  <0 


d, £ 

(25) 


d<Z 

(26) 

*0).  The 

(27) 
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P+  = -2sgn(/)7^(0,/l)^(0)/  f2 


+ 


sgn(/) 

f2  ' 


a 

x /K(cr,A)sy 

,  „  d^K^A) 

dcr 

Sy  da 

cos  [F(<7,  A)] 


+ 


)k/W(<X,A),l  dJK^~y l) 

- 1"  s  - 

dcr  do 


sin  [F(cr,  A)\ 


exp[F(cr,  /t)]c/cr 


(28a) 


P~  =2sgn(/)^(0, 2)^(0)/  f 


+ 


sgn(/)  ‘ 

f2  ' 


d 

y/K(<r,A)sy 

,  „  d^K(^A) 

do 

Sy  do 

sin  [F(cr,  2,)] 


d  yjK{(r,A)sx  dJK(^~A) 

- 1"  s  - 

da  dcr 


cos  [F(cr,  /l)] 


exp[F(cr,  A  )]dcr 


(28b) 


F(a,A) 


(29) 


Here,  sgn(/)  is  the  sign  function. 

The  nondimensional  barotropic  ( uE,vE ),  and  baroclinic  (AuE,AvE)  components  of  the 

Ekman  spiral  are  calculated  for  the  global  oceans  except  the  regions  near  the  equator  (5°S  - 
5°N)  using  (24a),  (24b),  (25a),  (25b)  with  the  monthly  mean  global  ocean  density  p  (kg  m"3) 
calculated  from  the  World  Ocean  Atlas  2009  temperature  and  salinity  data  (1°  x  1°  resolution) 
(http://www.nodc.noaa.gov/OC5/WOA09/pubwoaQ9.htmn  using  the  International 
Thermodynamic  Equation  of  Seawater  (http://www.teos- 1 0.org/pubs/TEOS- 1 0  Manual.pdf),  the 
computed  KPP  eddy  viscosity  data  K(cr,  X),  the  ocean  friction  velocity  (u*)  data  (1°  x  1°)  from 
http://iridl.ldeo.columbia.edU/SOURCES/.DASILVA/.SMD94/.climatologv/  (DaSilva  et  al. 
1994),  the  monthly  mixed  layer  depth  ( h )  data  from 

http://www.nodc.noaa.gov/OC5/WOA94/mix.html  (Monterey  and  Levitus  1997),  and  the  angle 
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of  surface  wind  (6)  data,  which  is  computed  from  the  monthly  zonal  wind  data  downloaded  from 
the  website:  http://iridUdeo.columbia.edU/SOURCES/.DASILVA/.SMD94/.climatology/.u3/  and 
the  monthly  meridional  wind  data  downloaded  from  the  website: 
http://iridl.ldeo.columbia.edU/SOURCES/.DASILVA/.SMD94/.climatology/.v3/. 

Fig.  5  shows  examples  of  dimensional  Ekman  spirals  (uE  =VEuE)  with  (solid  curve) 

and  without  (dashed  curve)  baroclinic  components.  The  upper  left  panels  (a),  and  (b)  show  the 
Ekman  spirals  at  Location-1  (11°N,  159°W)  (i.e.,  north  equatorial  Pacific)  with  large  baroclinic 
components.  The  upper  right  panels  (c),  and  (d)  indicate  the  Ekman  spirals  at  Location-2  (43°N, 
169°E)  (i.e.,  northwestern  Pacific  mid-latitude)  with  small  baroclinic  components.  Profiles  of  the 
horizontal  density  gradient  (dp  /  dx  ,dp  /  dy)  are  much  larger  at  Location- 1  (lower  left  panels) 
than  at  Location-2  (lower  right  panels). 

6.  Baroclinic  Effect 

The  baroclinic  portion  of  the  Ekman  spiral  can  be  effectively  determined  by  the  ratio  of 

the  vertical  integration  of  baroclinic  Ekman  component  over  the  Ekman  velocity, 

00 

J  J&4+A^da 

M=\  - •  (30) 

jyjul+vldcr 

0 

Horizontal  distribution  of  M  has  strong  seasonal  and  spatial  variability  with  large  M- values  (> 
0.2)  occurring  in  the  tropical  North  Pacific  Ocean,  tropical  Atlantic  Ocean  (10°N  -  25°N),  and 
eastern  Arabian  Sea  with  the  largest  value  of  0.9  in  the  central  tropical  North  Pacific  Ocean  near 
the  dateline,  and  with  small  M-values  (<  0.2)  occurring  in  the  southern  hemisphere  in  January, 
and  vice  versa  in  July  (Fig.  6).  Comparison  between  Fig.  6  and  Fig.  1  shows  negative  correlation 
between  2  and  M:  large  (small)  X  corresponds  to  small  (large)  M.  Such  negative  correlation  is 
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found  in  the  scatter  diagrams  of  (k,  M)  for  the  global  oceans  in  January  (Fig.  7a)  and  July  (Fig. 
7b)  with  linear  regression  equations, 

M  =  0.0239-0.0195/1  (R2  =0.561)  for  January,  (31) 

M  =  0.0305 -0.017U  (R2  =0.463)  for  July  .  (32) 

The  two  regression  equations  are  significant  on  the  level  of  0.0005  with  the  numbers  of  paired 
data  are  6945  in  January  (Fug.  7a)  and  6940  in  July  (Fig.  7b).  The  negative  correlation  between 
X  and  M  may  be  related  to  the  increase  of  the  KPP  eddy  viscosity  with  the  decrease  of  X 
especially  for  X  <  0. 

The  baroclinicity  parameter  is  identified  by  the  vertical  integration  of  the  magnitude  of 
horizontal  s  -  gradient  (crossing  the  mixed  layer)  scaled  by  f, 

oo 

b=-p\4s'+slda’  <33> 

J  0 

which  shows  evident  spatial  variability  and  weak  seasonal  variability  (Fig.  8).  Since  B  is 
inversely  proportional  to  f,  the  5-value  is  usually  large  (5  >  5)  in  low  latitudes  (20°S  -  20°N), 
and  small  ( B  <5)  in  middle  and  high  latitudes.  It  is  noted  that  the  scale  factor  of  f  for  the 
baroclinicity  parameter  (5)  [see  (33)]  (scaled  by  f)  is  only  used  for  searching  for  baroclinic 
Ekman  components  since  ( AuE,AvE )  are  inversely  proportional  to/  [see  (25),  (26),  (28a),  and 
(28b)]. 

The  baroclinic  effect  on  the  Ekman  spiral  is  evaluated  by  the  correlation  coefficient  (. R ) 
between  the  two  parameters,  M  and  b,  under  the  wind  (A>—1  and  combined  ( X  < -1) forcing 
regimes  in  the  northern  and  southern  hemisphere  (Table  1).  It  is  found  that  M  and  b  are 
positively  correlated  for  all  the  situations  with  large  R  (>  0.73)  for  both  hemispheres  and  months 
(January  and  July)  under  the  surface  wind  forcing  regime  (A>-1,  with  large  R  in  the  northern 
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hemisphere  in  January  (0.77)  and  in  the  southern  hemisphere  in  July  (0.62)  under  the  combined 
forcing  regime  (A  <  -1) ,  and  with  small  R  in  the  northern  hemisphere  in  July  (0.49)  and  in  the 
southern  hemisphere  in  January  (0.36). 

The  scatter  diagrams  of  ( b ,  M)  for  the  northern  hemisphere  (Fig.  9)  and  southern 
hemisphere  (Fig.  10)  also  show  the  similar  statistical  relationships  between  M  and  b  (M  increase 
as  B  increases).  Since  M  vanishes  as  b  vanishes,  (i.e.,  no  baroclinic  Ekman  components  when  the 
horizontal  density  gradient  equals  zero),  the  linear  regression  equation  between  M  and  b  is 
written  by 

M  =  cb ,  (34) 

where  the  regression  coefficient  c  is  obtained  using  the  least  square  error  method.  The 
regression  coefficient  c  is  always  positive  (Table  1).  It  has  largest  value  under  the  combined 
forcing  regime  in  January  for  the  northern  hemisphere  (0.0498),  and  in  July  for  the  southern 
hemisphere  (0.0370).  It  has  smallest  value  under  the  combined  forcing  regime  in  July  for  the 
northern  hemisphere  (0.00551)  (prevailing  wind  forcing  regime),  and  in  January  for  the  southern 
hemisphere  (0.00932)  (prevailing  wind  forcing  regime).  Thus,  the  baroclinic  effect  is  enhanced 
in  the  hemisphere  with  prevailing  combined  forcing  regime  and  weakened  in  the  hemisphere 
with  prevailing  wind  forcing  regime. 

7.  Eddy  Viscosity  due  to  Pure  Wind  Forcing 
7.1.  General  Description 

Earlier  studies  such  as  in  McWilliams  and  Huckle  (2006)  and  Song  and  Huang  (2011) 
assume  no  surface  buoyancy  flux  (B  =  0),  i.e.,  the  depth  ratio  X  =  0  [see  (17)],  the  depth 
dependent  non-dimensional  turbulent  velocity  scale  wx(o)  equals  1  [see  (20)  and  (21)].  Also,  the 
dimensional  form  of  the  Ekman  equation  is  used 
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ykxu. 


8z 


oz 


1  3  r  a  i 

7kx&M- 


with  the  surface  (dimensional)  boundary  condition, 

'da' 


K(  0) 


V  dz  yz=0 


=  t/J,+kx 


^W(0)S(0)A 


V 


/ 


and  the  lower  boundary  condition, 

u£  finite  as  z— »- oo, 


where 


if(0)  =  0.004w2/|/|. 


(35) 


(36a) 


(36b) 


(36c) 


With  the  monthly  mean  surface  wind  stress  data,  the  ocean  friction  velocity  ut  [using  (12b)]  and 


in  turn  the  surface  eddy  viscosity  K(0)  [using  (36c)]  are  calculated  except  the  equatorial  region 

5°S  -  5°N.  Fig.  1 1  clearly  shows  strong  horizontal  and  seasonal  variations  of  K(0)  .  In  January,  it 
has  large  values  (>  0.02  m2  s'1)  in  the  westem/central  tropical  North  Pacific  and  Atlantic  oceans 
(6°N  -  25°N),  medium  values  (0.01  -  0.02  m2  s'1)  in  the  mid-latitudes  associated  with  the 
Kuroshio  and  Gulf  Stream,  and  small  values  (<  0.01  m2  s'1)  in  rest  of  the  global  oceans. 
However,  in  July,  it  has  very  large  values  (>  0.03  m2  s'1)  in  the  western  Arabian  Sea  (related  to 
the  Southwest  monsoon),  large  values  (0.02  -  0.03  m2  s'1)  in  the  western  Bay  of  Bengal  and  the 
southern  tropical  Indian  and  Atlantic  oceans  (6°S  -  25°S),  medium  values  (0.01  -  0.02  m2  s'1)  in 
the  southern  tropical  Pacific  Ocean  (6°S  -  25°S),  mid-latitude  (35°S  -  45°S)  Indian  Ocean,  and 
small  values  (<  0.01  m2  s'1)  in  rest  of  global  oceans  in  July. 

The  eddy  viscosity  K(z )  has  three  different  types:  (a)  wind  and  depth  dependent  using 
the  KPP,  (b)  wind  dependent  and  depth  independent,  i.e.,  taking  surface  value  W(0)for  the 
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whole  water  column,  and  (c)  wind  and  depth  independent,  i.e.,  assigning  a  constant  value. 
Correspondingly,  the  solutions  for  the  three  types  of  eddy  viscosity  are  represented  by 


uf  +Auf,  vf  =v®+Avf,  i  =  1,2,3  (37) 

where  (u£\  v^°)  are  the  components  of  the  Ekman  velocity  in  barotropic  ocean  [i.e.,  the  Ekman 
solutions  when  (11)  is  satisfied];  and  (Am^Av^)  are  the  baroclinic  components  of  the  Ekman 
spiral.  The  baroclinic  effect  is  identified  by  the  root-mean  square  (RMS)  within  the  ocean  mixed 
layer  of  the  baroclinic  components  [( A ,  Av£}),  ( A ,  Av^2) ),  ( A ,  Avf ) )], 

R:  =  +[A^'0')]2}.  >-  =  1,2,3  (38) 

where  j  denotes  the  vertical  level;  and  J  is  the  total  number  of  the  vertical  levels  from  the 
surface  to  the  mixed  layer  depth. 

7.2.  Wind  and  Depth  Dependent  Eddy  Viscosity 
The  vertically  varying  eddy  viscosity  due  to  the  surface  wind  stress  is  given  by  (Song  and 
Huang  2011) 

k(z)  =  £(0X1  -  axz)  exp(a2z) ,  (39) 

where  (at,a2)  are  positive  constants.  Fitting  (37)  with  the  flow  in  the  /-plane  using  the  large- 

eddy  simulations  (Zikanov  et  al.  2003)  gives  the  following  semi-empirical  formula  (Song  and 
Huang  2011) 

K(z)  =  £(0)G(M-),  G(t)=(l - 64.0327/ exp(4.0073t)  ,  (40) 

u* 
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to  calculate  the  depth-dependent  eddy  viscosity  due  to  the  surface  wind  stress.  Here,  G{t )  is  the 
shape  function.  Fig.  12  shows  the  dependence  of  G{t )  versus  t,  where  t  is  the  non-dimensional 
depth,  t  =  \f\z/u+.  It  is  noted  that  ^(0)  is  inversely  proportional  to  the  magnitude  of  the 
Coriolis  parameter  ]/[.  For  the  same  ocean  friction  velocity  u *,  the  lower  the  latitude,  the  higher 
the  value  of^(0) . 

The  Ekman  velocity,  =\u^\v^]  are  the  approximate  analytical  solutions  of  (35)  by 

the  WKB  method  (see  Appendices  A  and  B)  with  the  eddy  viscosity  K(z )  given  by  (40).  The 
barotropic  components  are  given  by 

u\ ^  =  exp  [F(z)]  jc+ cos  [F(z)]  +  F_  sin  [F(z)]j,  (41a) 

vf  =exp[F(z)]{F+sin[F(z)]-F'cos[F(z)]},  (41b) 

where 


_  liTTf  k 

'  2  -#(o 


f  =  l'E<i  sgn(/K 

The  baroclinic  components  are  given  by 
A  uf  =  exp  [F(z)]  { A/+  cos  [F(z>]  +  \V~  sin  [F(z)]} 
+ 


+ 


+ 


[j,(0)-sy(0)]J 

2 

dK(£)  ^ 

4/  ii 

\f\K{()  \  S( 

[sx(0)  +  sy(0)]<j  1 

2 

dK(£)„ 

4/  ti 

\f\k(S)  \  B( 

[Sx(0)-Sym  r 

1  2 

I  dK(0 

4/  l\ 

||/|^(f)l  d( 

[sx(o)+sym  f 

1  2 

1  mo 

4/  L\ 

\\f\k(o  \  8( 

(42) 


cos[F(2)-F(0]lrap(-[F(r)-F(f)])</f 
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Av^  =  exp [F(z)] | AjF+  sin [F(z)]  - AtF  cos [F(z)]  j 


+ 


+ 


+ 


[5,(0)-5„(0)U 

2  \dk(0" 

4/  iv 

1/lAol  Sf 

[5,(0)  +  s,(0)U  I 

2  ax(Oc 

4/  J-V 

\f\k(S)\  d( 

[5,(0)  +  5,(0)]  J 

1  2  f  dk(C) 

4/  L\ 

\\f\k(C)\  K 

[Sx(0)-Sym  zr 

I  2  \dk(C) 

4/  L\ 

(|/|Aol  sc 

Here 


a  ,r=- 


pW  V<°>  M0>  ?exp[F(z)]fa^(z)  | 

V21/1  /  242 ml  4W)  i  *  1  (  )Jr 

i(°)  |expgz)]|s^z)iin[F(z)]^ 


2yl2f^\f\i  JI(z)  1  dz 


(44a) 


ffi0)  sx( 0)  ,  j,(0)  °f  exp[F(z)]  f dk(z) 


A  y-  =  r-^^^+ 


2|/|  /  272/J^I  JI(j)  1  5z 

^(Q)  f  exp[F(z)]|a^(z)  J 

27 l/^-l  I  5z  SinLF(z)Jj 


cos 


N] 


(44b) 


Fig.  13  shows  the  global  horizontal  distribution  and  zonal  mean  R\  (m  s'1)  in  January  and 
July.  In  January,  it  has  large  values  (>  0.2  m  s'1)  in  the  tropical  North  Pacific  and  Atlantic 
oceans  (6°N  -  20°N),  the  tropical  North  Indian  Ocean  (6°N  -  10°N),  the  tropical  South  Pacific 
and  Atlantic  oceans  (6°S  -  10°S),  medium  values  (0.02  -  0.2  m  s'1)  in  the  mid-latitudes 
associated  with  the  Kuroshio  and  Gulf  Stream,  south  tropical  Indian  Ocean  (6°S  -30°S)  and  small 
values  (<  0.01  m  s'1)  in  rest  of  the  global  oceans.  In  July,  it  has  large  values  (>  0.2  m  s'1)  in  the 
tropical  oceans  (15°S  -  15°N)  with  a  maximum  value  of  0.62  ms'1  in  the  Arabian  Sea.  It 
generally  decreases  with  the  increasing  latitude.  In  high  latitudes,  it  is  very  small  (<  0.01  m  s'1). 
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Zonal  mean  R\  decreases  with  latitude  near-exponentially  in  both  northern  and  southern 
hemispheres)  from  0.6  m  s"1  (0.4  m  s"1)  at  6°N  (6°S)  to  less  than  10~3  m  s"1  at  high  latitudes  near 
60°N  (60°S).  The  histograms  and  associated  probability  density  functions  (PDFs)  for  January  and 
July  (Fig.  14)  show  near-exponentially  decreasing  probability  with  R\  (log  scale  used  in  the 
vertical  axis).  The  probability  of  R\  larger  than  0.2  m  s"1  is  2.6%  [=1  -  P(R< 0.2  m  s"1)  =  1  - 
0.974)]  in  January  and  4.6%  (=  1  -  0.954)  in  July.  The  probability  of  R\  larger  than  0.1  m  s_1  is 
10.0%  in  January  and  1 1.5%  in  July.  The  95th  percentile  is  0.173  m  s'1  in  January  and  0.212  m  s"1 
in  July  (also  see  Tables  2  and  3). 

7.3.  Wind-dependent  and  Depth-independent  Eddy  Viscosity 
The  eddy  viscosity  is  given  by 


K(z)  =  K(  0). 

Substation  of  (45)  into  (42)  leads  to, 


F(z)  = 


a 


A 


w 


w 


2K(0)  _  1 0.008w* 

i/i  i  i/r 


(45) 


(46) 


where  (40)  is  used;  Dw  is  the  e-folding  decay  scale  of  the  Ekman  depth,  which  varies  with  the 
surface  wind  stress  through  /(0) ,  and  latitude.  For  the  same  ocean  friction  velocity  ( ut ).  the 
lower  the  latitude,  the  higher  the  value  of  Dw.  The  e-folding  scale  Dw  is  computed  from  the 
global  data  of  K(0)  using  (46).  It  also  has  evident  horizontal  and  weak  seasonal  variations  (Fig. 
15,  equatorial  region  5°S  -  5°N  not  computed):  large  values  (>  100  m)  generally  occur  in  sub¬ 
tropical  regions  (10°N  -  20°N,  10°S  -  20°S)  in  January  and  July  due  to  small  Coriolis  parameter. 
In  the  Northern  Hemisphere  middle  and  high  latitudes,  Dw  is  larger  in  January  (50-100  m  in  the 
Atlantic  Ocean  and  western  Pacific  Ocean)  than  in  July  (mostly  less  than  20  m).  In  the 
Southern  Hemisphere  middle  latitudes,  zonal  belts  with  medium  values  of  Dw  (20-50  m)  appear 
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at  35°S  -  50°S  in  the  three  southern  oceans  (Atlantic,  Indian,  and  Pacific)  in  January  but  does  not 
appear  in  the  Southern  Pacific  in  July. 

Substation  of  (46)  into  (41a),  (41b),  (43a),  (43b),  (44a)  and  (44b)  leads  to  the  barotropic 
components 


u(f2)  =  exp(— ) 


W 


D, 


V+  cos(— — )  +  V  sin(-^-) 
D  D 

^ w  ^ w 


(47a) 


vi2)  =  exp(— ) 


W 


D, 


V+  sin(-^-)  -  V  cos(— — ) 

D  D 

‘-'w 


and  the  baroclinic  components 


At/ 7  =  exp(^^) 


D, 


w 


A/+  cos(— — )  +  A/-  sin(— ) 


D, 


W 


D, 


W 


(47b) 


(48a) 


Av^2)  =  exp(— J-) 


D, 


W 


A/+  sin(— )  -  A/“  cos(— — ) 


D, 


W 


D, 


w 


(48b) 


Here  (V+,V  )  are  given  by  (42)  and 


A2V+  =  - 


A^(°)  AT>-_  Dwsx(0) 


f 


a2v~=- 


2/ 


(49) 


The  global  horizontal  distribution  of  R2  (m  s'1)  (Fig.  16)  is  similar  to  that  of  R\  (Fig.  13) 
in  both  January  and  July  with  smaller  values  in  the  Gulf  Stream  and  Kuroshio  regions. 
Latitudinal  decrease  of  zonal  mean  RMS  with  the  eddy  viscosity  is  also  comparable  in  both 
Northern  and  Southern  Hemispheres.  The  histograms  and  associated  PDFs  for  January  and  July 
(Fig.  17)  also  show  near-exponentially  decreasing  probability  with  R2.  For  large  Rj  values  (>  0.2 
m  s'1),  the  probability  is  zero  in  January  and  0.6%  (=  1  -  0.994)  in  July.  The  probability  of  R2 
larger  than  0.1  m  s'1  is  0.28%  in  January  and  4.7%  in  July.  The  95th  percentile  is  0.089  m  s'1  in 
January  and  0. 1 1  m  s'1  in  July. 
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7.4.  Wind  and  Depth  Independent  Eddy  Viscosity 


For  the  wind  and  depth  independent  eddy  viscosity, 

K(z)  =  K0=  const .  (50) 


Substitution  of  (50)  into  (46)  gives 


DW=D  = 


(51) 


which  is  the  classical  Ekman  depth.  After  replacing  Dw  by  D,  (47a)  and  (47b)  lead  to  the 
barotropic  components  (u^  ,v^])  and  (48a)  and  (48b)  lead  to  the  baroclinic  components 


(  Am^3)  ,  Av^3)  )  of  the  Ekman  velocity.  Here,  the  constant  eddy  viscosity  K{)  is  taken  as  0.054 
m2  s'1  (Price  et  al.  1987). 

The  global  horizontal  distribution  of  R3  (m  s'1)  (Fig.  18)  is  similar  to  that  of  R]  (Fig.  13) 
in  January  and  July  with  latitudinal  decrease  of  zonal  mean  RMS  in  the  Northern  and  Southern 
Hemispheres.  The  histograms  and  associated  PDFs  for  January  and  July  (Fig.  19)  also  show 
near-exponentially  decreasing  probability  with  R}.  For  large  R}  values  (>  0.2  m  s'1),  the 
probability  is  zero  in  January  and  July.  The  probability  of  R}  larger  than  0.1  m  s'1  is  0.1%  in 
January  and  1.2%  in  July.  The  95th  percentile  is  0.04  m  s'1  in  January  and  0.058  m  s'1  in  July. 

8.  Conclusions 


Analytical  solution  of  the  Ekman  spiral  in  real  oceans  is  obtained  with  vertical 
geostrophic  and  ageostrophic  shears  linking  to  turbulent  stress  in  upper  oceans,  under  surface 
wind  and  buoyancy  forcing  using  the  KPP  eddy  viscosity.  The  Ekman  spiral  contains  barotropic 
and  baroclinic  components.  The  barotropic  component  is  similar  to  the  classical  Ekman  spiral. 
The  baroclinic  component  is  caused  by  horizontally  inhomogeneous  density.  The  baroclinic 
component  vanishes  as  the  horizontal  density  gradient  vanishes. 
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Climatological  monthly  data  of  global  ocean  mixed  layer  depth,  Monin-Obukhov  length 
scale,  friction  velocity,  surface  winds,  and  density  profiles  are  used  to  calculate  the  depth  ratio 
(2),  KPP  eddy  viscosity,  the  barotropic  and  baroclinic  Ekman  velocities,  the  baroclinicity 
parameter  (5),  and  the  proportion  of  the  baroclinic  Ekman  component  (M).  Large  baroclinic 
proportion  is  usually  associated  with  the  prevailing  combined  forcing  regime  such  as  in  the 
northern  (southern)  hemisphere  in  January  (July). 

Statistical  analysis  on  the  calculated  global  (2,  M,  B )  values  shows  significant  negative 
correlation  between  2  and  M:  large  (small)  X  corresponds  to  small  (large)  M,  and  significant 
positive  correlation  between  B  and  M:  large  (small)  B  corresponds  to  large  (small)  M.  The 
negative  correlation  between  A  and  M  may  be  related  to  the  increase  of  the  KPP  eddy  viscosity 
with  the  decrease  of  X  especially  for  A  <  0.  The  positive  correlation  coefficient  between  B  and  M 
varies  with  the  prevailing  wind  and  combined  forcing  regimes.  The  baroclinic  effect  is  enhanced 
in  the  hemisphere  with  prevailing  combined  forcing  regime  and  weakened  in  the  hemisphere 
with  prevailing  wind  forcing  regime. 

For  pure  wind  forcing  (i.e.,  zero  surface  buoyancy  flux),  three  types  of  eddy  viscosity 
from  existing  parameterization  [wind  and  depth  dependent,  wind-dependent  and  depth- 
independent,  and  wind  and  depth  independent  (i.e.,  constant  eddy  viscosity)]  as  well  as  the 
vertical  root-mean  square  of  the  baroclinic  component  within  the  ocean  mixed  layer  (R)  of  the 
analytical  Ekman  spiral  are  used  to  investigate  the  baroclinic  effect.  It  enhances  with  the 
decreasing  latitude  and  usually  very  evident  (>  0.2  m  s"1)  in  the  tropical  oceans  in  January  and 
July,  and  extremely  large  value  of  0.62  m  s"1  in  the  Arabian  Sea  in  July  for  all  the  cases  (three 
types  of  eddy  viscosity).  In  middle  and  high  latitudes  (especially  in  the  South  Hemisphere),  it  is 
generally  very  small  (i.e.,  the  classical  Ekman  spiral  applies)  except  in  the  Gulf  Stream  and 
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Kuroshio  regions  in  January.  These  results  are  consistent  with  the  earlier  observational  studies 
such  as  conducted  in  the  Drake  Passage  (Lenn  and  Chereskin  2009). 

The  near-exponentially  decreasing  probability  with  the  vertical  root-mean  square  of  the 
baroclinic  component  is  obtained  from  the  histograms.  The  statistical  characteristics  show  that 
the  baroclinic  components  for  the  three  types  of  eddy  viscosity  under  pure  wind  forcing  are  all 
comparable  in  January  and  July.  Near-exponentially  decreasing  probability  with  R2 ,  or  R2  is 
found.  The  probability  of  R\  larger  than  0.2  m  s"1  is  2.6%  [=1  -  P(R<  0.2  m  s"1)  =  1  -  0.974)]  in 
January  and  4.6%  (=  1  -  0.954)  in  July.  The  probability  of  R\  larger  than  0.1  m  s"1  is  10.0%  in 
January  and  11.5%  in  July.  The  95th  percentile  is  0.173  m  s"1  in  January  and  0.212  m  s"1  in  July 
(also  see  Tables  2  and  3). 

Finally,  it  is  noted  that  the  monthly  mean  density  fields  from  the  WOA-2009  temperature 
and  salinity  data  cannot  represent  density  fronts  associated  with  submesoscale  processes.  The 
computation  here  is  only  to  show  the  importance  of  horizontally  inhomogeneous  density  on  the 
Ekman  spiral.  Further  computation  is  needed  to  verify  the  good  approximate/analytical  solutions 
again  high  horizontal-resolution  wind  and  density  data  if  they  will  be  available. 

Appendix-A  General  Solutions  of  Eq(22) 

Let  the  Ekman  currents  ( uE ,  vE)  be  represented  by  a  complex  variable  y/, 


(Al) 


Substitution  of  (Al)  into  (22)  leads  to 


(A2) 


where 


(A3) 
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The  function  $  represents  the  baroclinic  effect. 


The  second-order  differential  equation  (A2)  needs  two  boundary  conditions.  The  surface 
boundary  condition  (14)  becomes 


dy/ 


do  a-°  kK(0)Vl 


U*  (cos  0  +  i sin  0) +  „  -  -  |>v (0)  - isx (0)] . 


(A4) 


The  lower  boundary  condition  of  equation  (A2)  is  given  by 

|^|  finite  as  cr— » oo. 


(A5) 


to  guarantee  a  physically  meaningful  solution,  i.e.,  \y/\  cannot  be  infinity  as  o  — >co. 

Eq.  (A2)  is  a  linear  inhomogeneous  ordinary  differential  equation  with  the  depth-varying 
coefficient  K(z).  Following  Berger  and  Grisogono  (1998),  studying  the  Ekman  atmospheric 
boundary  layer,  an  approximate  solution  to  the  inhomogeneous  problem  (A2)  can  be  found 
with  the  variation  of  parameters  technique,  provided  that  an  approximate  solution  of  the 
homogeneous  problem  of  (A2), 


_  _  d  V  dK  dy/  .  . 

K — V  + - —  -i  ft¥  =  0 , 

do 2  do  do 


(A  6) 


exists.  If  two  independent  approximate  solutions  to  homogeneous  problem  of  Eq.  (A6)  are  given 
by  y/\  (o)  and  i//2  (o) ,  the  general  solution  of  Eq.  (A2)  is  given  by 


¥  =  <\¥ i  (o')  +  ci¥i  (cr)  +  c)  (o)y/,  ( o)  +  c2  (o)y/2  (cr)  , 


(A7) 


where 


^(o)  =  -J 


yA 


K(OW,  (Qdy  2  (O I  d£-v  2  (Od¥i  (C)  /  dQ 


dC, 


(A8) 


ou 


y/x(j) 


i  K(0[¥x (Qdy 2 (O / dC- ~ ¥2 (QdVx (O / d£J 


d£ . 


(A9) 
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Appendix-B  The  WKB  Method  for  Solving  Eq.(A6) 

The  WKB  method  can  be  used  to  obtain  a  good  approximate  solution  of  Eq.  (A6)  if  the 


vertical  variation  of  K(a)  is  slower  than  that  of  y/((j)  (Grisogono  1995), 


iff  oc  exp 


(S0+Sl£  +  S2£2+...) 


where  s  is  a  presumably  small  parameter.  Substitution  of  (Bl)  into  (A6)  leads  to  < 
equations  in  terms  of  powers  of  s .  If  K(g)  does  not  vary  too  quickly  with  depth,  we  have 


L. - n  =  0,  1,  2,... 

iwl 


Solving  for  the  first  two  terms  So  and  S\  yields 


>■ 


5,  =  —  In 


^(0) 

K(a) 


Thus,  the  two  approximate  solutions  of  the  homogeneous  equation  (A4)  are 
¥i (o’)  =  A<r)  exp[(l  +  i)F(cr)\,  i//2 (cr)  =  A(cr)  exp[-(l  +  i)F(cr)] , 


where 


Act)  = 


Kj  0) 


nl/4 


,  Fip’)  =  ~J~~j 


2  J0 


Substitution  of  (B5)  and  (B6)  into  (A8)  and  (A9)  gives 


{&)  =  - 


(1-0 


2^2|/0|X(0)  { 


jA(0<HOexp[-O  +  mO]d£, 


c2(cr)  =  - 


(1-0 


2^21/01^(0)  i 


oo 

jA(0fi(Oexp[(l+i)F(0\d£. 


(Bl) 

set  of 

(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 
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Substitution  of  (B5)  into  (A7)  gives 

¥  =  \Pi+  G  {o)\A{a)  exP[(l  +  O^(o-)]  +  [c2  +  c2  (cr)]A(cr)  exp[-(l  +  i)F(a)\ .  (B9) 

It  is  noted  that  F(a)  <  0  leads  to 

y/2  — »  oo  as  <j  — >  oo .  (BIO) 

This  leads  to 


c2  =  0  (Bll) 

Substitution  of  (B9)  into  the  surface  boundary  condition  (A4)  gives 


G 


c2(0)  +  (l-0 


w.  (cos#-!-/ sin#) 


K(0)  h 

|2/0|  / 


[sy(o)-isxm\- 


=  c2(0)  +  V+  —iV~  + 


\K( 0)  h 

'M7 


(°)  -  Sy  (0)]  +  i  [_sx  (0)  +  sy  (0)]} . 


(Bll) 
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Table  1 .  January  and  July  correlation  coefficients  and  number  of  paired  data  between  (M,  B )  for 
northern  and  southern  hemispheres  under  wind  and  combined  forcing  regimes. 


Location 

Statistics 

January 
(A  >  -1) 

January 

(A<-1) 

July 
(A  >  -1) 

July 

(A<-1) 

Northern 

Hemisphere 

Number 
of  Paired 
Data 

1181 

1225 

2120 

143 

R 

0.753 

0.773 

0.797 

0.487 

c 

0.00939 

0.0498 

0.0101 

0.00551 

Southern 

Hemisphere 

Number 
of  Paired 
Data 

4015 

220 

2573 

1660 

R 

0.794 

0.358 

0.733 

0.624 

c 

0.0142 

0.00932 

0.0134 

0.037 

Table  2.  Statistical  characteristics  of  the  VRMA  (within  the  ocean  mixed  layer)  of  the  baroclinic 
components  over  the  global  oceans  for  the  three  types  of  eddy  viscosity  under  zero  surface 
buoyancy  flux  in  January. 


RMS 

Qo.5  (m/s) 

Q0.95  (m/s) 

P(R<  0.02  m/s) 

P(R<0.1  m/s) 

P(R<0.2  m/s) 

Ri 

0.0113 

0.173 

0.669 

0.900 

0.974 

Ri 

0.00401 

0.0888 

0.800 

0.972 

1.000 

Ri 

0.00107 

0.0399 

0.901 

0.999 

1.000 

Table  3.  Statistical  characteristics  of  the  VRMA  (within  the  ocean  mixed  layer)  of  the  baroclinic 
components  over  the  global  oceans  for  the  three  types  of  eddy  viscosity  under  zero  surface 
buoyancy  flux  in  July. 


RMS 

Qo.5  (m/s) 

Q0.95  (m/s) 

P(R<  0.02  m/s) 

P(R<0.1  m/s) 

P(R<0.2  m/s) 

Ri 

0.00929 

0.212 

0.694 

0.875 

0.954 

Ri 

0.00306 

0.110 

0.795 

0.953 

0.994 

R3 

0.000769 

0.0575 

0.883 

0.988 

1.000 
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Figure  Captions 


Fig.  1.  Monthly  depth  ratio  (X):  (a)  January,  and  (b)  July.  It  is  noted  that  only  two  regimes  are 
evident  with  the  monthly  mean  data:  wind-forcing  and  combined  forcing  regimes.  Here,  the 
black  contours  are  referred  as  X  =  -1. 

Fig.  2.  The  shape  function  G(o). 

Fig.  3.  Dependence  of  (a)  (/)  and  (b)  K  on  a  and  X  for  1  >  cr  >  0 . 

Fig.  4.  Monthly  Ekman  velocity  scale  (Ve):  (a)  January,  and  (b)  July. 

Fig.  5.  Examples  of  Ekman  spirals  (a)  uE,  (b)  vE  at  Location-1  (11°N,  159°W),  (c)  uE,  and 
(d)  vE  at  Location-2  (43°N,  169°E)  (upper  panels)  as  well  as  corresponding  horizontal  density 

gradients  (lower  panels).  It  is  noted  that  the  horizontal  density  gradients  are  much  stronger  at 
Location- 1  than  Location-2. 

Fig.  6.  Monthly  horizontal  distribution  of  M:  (a)  January,  and  (b)  July. 

Fig.  7.  Scatter  diagrams  of  (2,  M)  with  linear  regression:  (a)  January,  and  (b)  July.  It  is  noted  that 
the  negative  correlation  between  (2,  M)  is  significant  on  the  level  of  0.0005. 

Fig.  8.  Monthly  horizontal  distribution  and  zonal  mean  of  B :  (a)  January,  and  (b)  July. 

Fig.  9.  Scatter  diagrams  of  (5,  M)  with  linear  regression  for  the  northern  hemisphere:  (a) 
January,  wind  forcing  regime  (A  >  -1,  (b)  January,  combined  forcing  regime  (A  <  -1),  (c)  July, 
wind  forcing  regime  (A  >  -1,  and  (d)  July,  combined  forcing  regime  (A  <  -1).  It  is  noted  that  the 
positive  correlation  between  (2,  M)  is  significant  on  the  level  of  0.0005. 

Fig.  10.  Scatter  diagrams  of  (5,  M)  with  linear  regression  for  the  southern  hemisphere:  (a) 
January,  wind  forcing  regime  (A  >  -1,  (b)  January,  combined  forcing  regime  (A  <  -1),  (c)  July, 
wind  forcing  regime  (A  >  -1,  and  (d)  July,  combined  forcing  regime  (A  <  -1).  It  is  noted  that  the 
positive  correlation  between  (2,  M)  is  significant  on  the  level  of  0.0005. 

Fig.  11.  Horizontal  distribution  of  eddy  viscosity  at  the  ocean  surface  A^(O)  calculated  from 
monthly  mean  surface  wind  stress  data  (downloaded  from  the  websites: 
http://iridl.ldeo.columbia.edU/SOURCES/.DASILYA/.SMD94/.climatology/.u3/  and 

http://iridl.ldeo.columbia.edU/SOURCES/.DASILVA/.SMD94/.climatology/.v3/ 

using  (20)  and  (21):  (a)  January,  and  (b)  July. 

Fig.  12.  Vertical  structure  function  G(t)  versus  t. 
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Fig.  13.  Horizontal  distribution  and  zonal  mean  of  vertical  root-mean  square  of  baroclinic 
components  of  the  Ekman  spiral  in  neutral  ocean,  R\  (m  s"1),  inside  the  mixed  layer  with  wind- 
dependent  and  depth-dependent  eddy  viscosity  K{z)\  (a),  January,  and  (b)  July. 

Fig.  14.  Histogram  and  probability  density  function  of  R\  (ms"1)  with  wind-dependent  and  depth- 
dependent  eddy  viscosity  K(z )  :  (a)  January,  and  (b)  July. 

Fig.  15.  Horizontal  distribution  of  e-folding  depth  ( Dw )  of  the  Ekman  layer  with  wind-dependent 
and  depth-independent  eddy  viscosity:  (a)  January,  and  (b)  July. 

Fig.  16.  Horizontal  distribution  and  zonal  mean  of  vertical  root-mean  square  of  baroclinic 
components  of  the  Ekman  spiral  in  neutral  ocean,  R2  (m  s"1),  inside  the  mixed  layer  with  wind- 

dependent  and  depth-independent  eddy  viscosity  A"(0)  :  (a),  January,  and  (b)  July. 

Fig.  17.  Histogram  and  probability  density  function  of  R2  (m  s"1)  inside  the  mixed  layer  with 
wind-dependent  and  depth-independent  eddy  viscosity  7^(0) :  (a),  January,  and  (b)  July. 

Fig.  18.  Horizontal  distribution  and  zonal  mean  of  vertical  root-mean  square  of  baroclinic 
components  of  the  Ekman  spiral  in  neutral  ocean,  R2  (m  s”1),  inside  the  mixed  layer  with  a 
constant  eddy  viscosity  (0.054  m2  s"1):  (a),  January,  and  (b)  July. 

Fig.  19.  Histogram  and  probability  density  function  of  R2  (m  s"1)  inside  the  mixed  layer  with  a 
constant  eddy  viscosity  (0.054  m2  s"1):  (a),  January,  and  (b)  July. 
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(a):  January 


Fig.  1.  Monthly  depth  ratio  (X):  (a)  January,  and  (b)  July.  It  is  noted  that  only  two  regimes  are 
evident  with  the  monthly  mean  data:  wind-forcing  and  combined  forcing  regimes.  Here,  the 
black  contours  are  referred  as  X  =  -1. 
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Fig.  2.  The  shape  function  G(o). 
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}  4  1  4 

A  xIO  A  xIO 

Fig.  3.  Dependence  of  (a)  </>  and  (b)  K  on  a  and  A,  for  1  >  cr  >  0 . 


38 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 


Latitude  Latitude 


(a):  January 


30 


(b):  July 
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Fig.  4.  Monthly  Ekman  velocity  scale  ( Ve ):  (a)  January,  and  (b)  July. 
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Location  1:  Lat=  11  Lon=  201  \=  -0.598 


-0.6  -0.4  -0.2  0  0.2  -0.4  -0.2  0  0.2  0.4 

uE  (m/s)  vE  (m/s) 


Location  2:  Lat=  43  Lon=  169  X=  -0.809 


-0.4  -0.2  0  0.2  0.4  -0.4  -0.2  0  0.2  0.4 

uE  (m/s)  vE  (m/s) 


Fig.  5.  Examples  of  Ekman  spirals  (a)  uE ,  (b)  vE  at  Location-1  (11°N,  159°W),  (c)  uE9  and 
(d)  vE  at  Location-2  (43°N,  169°E)  (upper  panels)  as  well  as  corresponding  horizontal  density 

gradients  (lower  panels).  It  is  noted  that  the  horizontal  density  gradients  are  much  stronger  at 
Location- 1  than  Location-2. 
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Latitude  Latitude 


(a):  January 


30  60  90  120  150  180  210  240  270  300  330 

Longitude 

^ _  ,  . 

0.02  0.04  0.06  0.08  0.1  0.12  0.14  0.16  0.18  0.2  0.22  0.24  0.26  0.28 

M 

Fig.  6.  Monthly  horizontal  distribution  of  M:  (a)  January,  and  (b)  July. 


(b):  July 


41 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 


M=0.0239-0.0195X;  R2=0.561  (6945  points) 

0.25 

0.2 

0.15 

0.1 

0.05 

0 

M=0.0305-0.0171X;  R2=0.463  (6940  points) 
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Fig.  7.  Scatter  diagrams  of  (2,  M)  with  linear  regression:  (a)  January,  and  (b)  July.  It  is  noted  that 
the  negative  correlation  between  (2,  M)  is  significant  on  the  level  of  0.0005. 
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Latitude  Latitude 


(a):  January 


30  60  90  120  150  180  210  240  270  300  330 

Longitude 


0.02  0.05  0.1  0.2  0.5  1  2  5  10 

B 

Fig.  8.  Monthly  horizontal  distribution  and  zonal  mean  of  B:  (a)  January,  and  (b)  July. 
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Fig.  9.  Scatter  diagrams  of  (5,  M)  with  linear  regression  for  the  northern  hemisphere:  (a) 
January,  wind  forcing  regime  (A  >  -1,  (b)  January,  combined  forcing  regime  (A  <  -1),  (c)  July, 
wind  forcing  regime  (A  >  -1,  and  (d)  July,  combined  forcing  regime  (A  <  -1).  It  is  noted  that  the 
positive  correlation  between  (A,  M)  is  significant  on  the  level  of  0.0005. 
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Fig.  10.  Scatter  diagrams  of  ( B ,  M)  with  linear  regression  for  the  southern  hemisphere:  (a) 
January,  wind  forcing  regime  (A  >  -1,  (b)  January,  combined  forcing  regime  (A  <  -1),  (c)  July, 
wind  forcing  regime  (A  >  -1,  and  (d)  July,  combined  forcing  regime  (A  <  -1).  It  is  noted  that  the 
positive  correlation  between  (A,  M)  is  significant  on  the  level  of  0.0005. 
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Fig.  11.  Horizontal  distribution  of  eddy  viscosity  at  the  ocean  surface  K( 0)  calculated  from 
monthly  mean  surface  wind  stress  data  (downloaded  from  the  websites: 
http://iridl.ldeo.columbia.edU/SOURCES/.DASILVA/.SMD94/.climatologv/.u3/  and 

http://iridl.ldeo.columbia.edU/SOURCES/.DASILVA/.SMD94/.climatologv/.v3/ 

using  (20)  and  (21):  (a)  January,  and  (b)  July. 
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Fig.  13.  Horizontal  distribution  and  zonal  mean  of  vertical  root-mean  square  of  baroclinic 
components  of  the  Ekman  spiral  in  neutral  ocean,  R\  (m  s"1),  inside  the  mixed  layer  with  wind- 
dependent  and  depth-dependent  eddy  viscosity  K{z)\  (a),  January,  and  (b)  July. 
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Fig.  14.  Histogram  and  probability  density  function  of  R\  (m  s’1)  with  wind-dependent  and  depth- 
dependent  eddy  viscosity  K(z) :  (a)  January,  and  (b)  July. 
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Fig.  15.  Horizontal  distribution  of  e-folding  depth  ( Dw )  (unit:  m)  of  the  Ekman  layer  with  wind- 
dependent  and  depth-independent  eddy  viscosity:  (a)  January,  and  (b)  July. 
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Fig.  16.  Horizontal  distribution  and  zonal  mean  of  vertical  root-mean  square  of  baroclinic 
components  of  the  Ekman  spiral  in  neutral  ocean,  R2  (m  s"1),  inside  the  mixed  layer  with  wind- 

dependent  and  depth-independent  eddy  viscosity  A^(0)  :  (a),  January,  and  (b)  July. 
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Fig.  17.  Histogram  and  probability  density  function  of  7?2  (m  s"1)  inside  the  mixed  layer  with 
wind-dependent  and  depth-independent  eddy  viscosity  £(0) :  (a),  January,  and  (b)  July. 
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Fig.  18.  Horizontal  distribution  and  zonal  mean  of  vertical  root-mean  square  of  baroclinic 
components  of  the  Ekman  spiral  in  neutral  ocean,  f?3  (m  s'1),  inside  the  mixed  layer  with  a 
constant  eddy  viscosity  (0.054  m2  s'1):  (a),  January,  and  (b)  July. 
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Fig.  19.  Histogram  and  probability  density  function  of  7?  3  (m  s'1)  inside  the  mixed  layer  with  a 
constant  eddy  viscosity  (0.054  m2  s'1):  (a),  January,  and  (b)  July. 
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